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ABSTRACT 

We present a sample of 23 spectroscopically confirmed Type la supernovae that were discovered in 
the background of galaxy clusters targeted by ROTSE-IIIb and use up to 18 of these to determine 
the local (z — 0.05) volumetric rate. Since our survey is flux limited and thus biased against fainter 
objects, the pseudo- absolute magnitude distribution (pAMD) of SNe la in a given volume is an 
important concern, especially the relative frequency of high to low-luminosity SNe la. We find that 
the pAMD derived from the volume limited Lick Observatory Supernova Search (LOSS) sample is 
incompatible with the distribution of SNe la in a volume limited (z < 0.12) sub sample of the SDSS- 
II. The LOSS sample requires far more low-luminosity SNe la than the SDSS-II can accommodate. 
Even though LOSS and SDSS-II have sampled different SNe la populations, their volumetric rates 
are surprisingly similar. Using the same model pAMD adopted in the SDSS-II SNe la rate calculation 
and excluding two high-luminosity SNe la from our sample, we derive a rate that is marginally higher 
than previous low-redshift determinations. With our full sample and the LOSS pAMD our rate is 
more than double the canonical value. We also find that 5 of our 18 SNe la are hosted by very 
low- luminosity (Mb > —16) galaxies, whereas only 1 out 79 nearby SDSS-II SNe la have such faint 
hosts. It is possible that previous works have under-counted either low luminosity SNe la, SNe la 
in low luminosity hosts, or peculiar SNe la (sometimes explicitly), and the total SNe la rate may be 
higher than the canonical value. 
Subject headings: supernovae: general 



1. INTRODUCTION 

Calculating the rate of Type la supernovae (SNe la) 
and its dependence on redshift illuminates an important 
contributor to the metal enrichment history of the uni- 
verse, the production rate of a specific cross section of 
stellar systems, and, when linked to star formation his- 
tory, it can reveal the nature and relative fractions of 
the progenitor systems. Theory and recent empirical ev- 
idence argues that the explosions are derived from de- 
generate, hydrogen depleted stars, namely white dwarfs 
(jHovle fc FowIerl[l96fl iNugent et al.ll20Tlh . The nature 
of the companion star which must be present to feed 
mass onto these otherwise stable constructions, however, 
is more controversial. 

It could very well be that there are two or more path- 
ways leading to what we observe as SNe la including 
events from single degenerate progenito rs (e.g. a white 
dwar f accreting from a red giant star; iWhelan fc Ibenl 
I1973D a nd double degen erate progenitors (i.e. two white 
dwarfs; iWebbinkl Il98l . This is a potential obstacle 
for the further use of SNe la as cosmological probes if 
the local population used to calibrate the peak magni- 
tude to light curve width relation differs systematically 
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from the distant, cosmologically significant population 
(i.e. if high redshift events obey a different relation, cos- 
mology studies may be biased; IDominguez et all 120011: 
i Mannucci et all 120061: ISullivan et al.l 120061 : iHowell et al l 
l2007t iQuimbv et al.ll2007a| ). 

The SNe la rate is one tool for resolving this issue. 
Locally, the quantity of progenitor systems demanded 
by the SNe la ra te may be compared to the actual sup- 
ply available Ce.g.|Nelemans et alJ [2005l: iKilic e t~aT1l2012l : 
2012h . At greater distances, the delay 
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time distribution (DTD), the distribution of SNe la pro- 
genitor systems lifespans from birth to explosion, can be 
compared to th e distributions expected fo rm various pro- 
genitor models (jYungelson fc Lividl2000l ) . The DTD can 
be recovered using the SNe la rate and st ar-formation 
histo ries of a targeted population of hosts (jMaoz et al.l 
2011), or by comparing SNe la over a r ange of redshifts 
to the cosmic star forma tion history (e.g. IMannucci et al.l 
120061 : iGraur et al]|2011[ ). The results are consistent with 
a power-la w distribution that favor s short delays (for a 
review, see iMaoz fc M annucci 20ll|). 

Indeed, SNe la appear to occur more frequently in late 
type gal axies with active star formation than in e llipti- 
cals (e.g. IMannucci et al.ll2005l: ISullivan et al.ll2006D . Re- 
cently, the Lick Observatory Supernova Search (LOSS) 
reported that the rate of SNe la per unit luminosity 
of the hosts is higher for low l uminosity galaxies than 
luminous ones (|Li et al.l l2011al ) . This so-called "rate- 
size" relation may be conne c ted to a metallicity effect 
(|Kobavashi fc Nomotol 12001 iKistler et all [20Tlft . since 
lower luminosity hosts tend to have lower global metallic- 
ities (jTremonti et al.l [2004) . White dwarfs that form in 
low metallicity environments may have a mass distribu- 
tion that is biased to more massive objects as compared 
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to solar metallicity populations (|Umeda et al.lll999l) , and 
this in turn may lead to a greater supply of SNe la pro- 
genitors. 

Measurement of the SNe la rate was originally per- 
formed through searches targeting specific galaxies, and 
thus the rates derived were in units of th e num- 
ber of SNe la found per time per galaxy (|Zwickvl 
1938). This evolved into a rate per galaxy luminosity- 
SNuB or SNuK, which give the number of SNe la per 
10 10 L p, (in the B or K-band s, respectively), per cen- 
tury. ICappellaro et al.l (|1999l ) combined targeted pho- 
tographic and visual surveys and found a SNe la rate 
of SNuB= 0.2 ± 0.06 averaged over all galaxy types. 
LOSS finds roughly compatible SNe la SNuB rates in 
specific host types if th ey do not include the rate-size 
relation (|Li et al.ll2011a| ). They a lso use an adopted K- 
band galaxy luminosity function ([Kochanek et al.l 1200 ll ) 
to convert their measurements into volumetric rates. 

Methods to dis cover SNe la without preference for 
host galaxies (e.g. lHamuv et al] 119931 : IPerlmutter et al.l 
1995), have been widely adopted as the availability of 
wide field cameras has increased (e.g. SNLS, ESSENCE, 
SDSS-II, PTF, PanSTARRS, SkyMapper). SNe la rates 
from such surveys are most readily reported in volumet- 
ric units, SNe IaMpc _3 yr _1 . Included in this survey 
category is the SDSS-II, which was designed to bridge 
the gap between the local SNe la and the high red- 
shift events used as cosmological probes (Fricman et al. 
2008). Over a three year period, the SDSS-II spectro- 
scopically con firmed several hundr ed SNe la in the range 
<; z <; 0.4. iDildav et all (|2010D use a sub sample of 
(normal) SNe la discoveries to derive the z < 0.3 rate, 
which is seen to increase roughly as a power-law with 
redshift, 1Z oc (1 + z) 2 . SNe la rate studies have been 
conducted at even higher redshifts, including work done 
with the Supernova Leg acy S urvey (SNLS; Neill et aTl 
l200i IPerrett et al1l20L¥. HST (IDahlen et al.ll2008lh and 
Subaru (|Graur et al.ll201lh . The measured SNe la rates 
increase to z ~ 1, and then level off. For a recent com- 
pil ation of volumetric SNe la rates from various surveys, 
see IGraur et al.l lj20H . 

In this paper, we determine the volumetric SNe la 
rate and contributions to this from dwarf and giant host 
galaxies based on discoveries from the ROTSE-IIIb tele- 
scope. Section [2] discusses how our sample was selected. 
We discuss our unfiltered magnitude system and how this 
may relate to SNe la pseudo-absolute magnitude distri- 
butions (pAMDs) constructed from filtered observations 
in Sj3] We also discuss the pAMDs derived from the LOSS 
and SDSS-II in this section. For the later, we deter- 
mine an empirical pAMD based on a volume limited sub- 
sample of the SDSS-II after discussing the pAMD model 
that was actually used in the SDSS-II SNe la rate calcu- 
lation (and which we adopt for comparison to the SDSS- 
II). In S]Hwe study our ability to detect point sources of 
various brightness as a function of seeing and limiting 
magnitudes, and we use this in §5] to determine our over- 
all survey efficiency-the probability of discovering SNe la 
drawn from a given population-as a function of distance. 
Next, we measure multi-band photometry for the host 
galaxies of our supernovae and calculate rest frame ab- 
solute magnitudes in ^ SNe la rates are calculated in 
and final conclusions are offered in £[8] Throughout 
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Fig. 1.— Winter-Spring search fields for the ROTSE-IIIb sur- 
veys. The color scale indicates the total number of nights each 
field was surveyed (i.e. reference epochs are not included). The 
diamonds mark the locations of background SNe la. 

this paper, we assume a flat, Ho = 71 kms -1 Mpc -1 , 
f2 m = 0.27 cosmology. 

2. SAMPLE 

Our sample is drawn from the background population 
of supernovae discovered by ROTSE-III b in the course 
of the Texas Supernova Search (TSS; IQuimbvl 120061) 
and the ROTSE S upernova Verification Project (RSVP; 
lYuan et al.l [2007al) . Details of the full ROTSE-IIIb su- 
pernova sample will be presented in a forthcoming paper. 
Here we summarize the key characteristics of the sample 
selection process. 

The TSS began survey operations in November of 2004 
and was succeeded by RSVP in early 2008. Our survey 
instrument, ROTSE-IIIb, has a 0.45 m aperture and a 
1.85 x 1.85 degree field of view. The search continues to 
operate and has been expanded to the other ROTSE-III 
telescopes in Australia, Namibia, and Turkey, but for this 
work we consider only discoveries made by ROTSE-IIIb 
in Texas prior to February 2009. 

We typically targeted galaxy clusters such as Virgo, 
Ursa Major, Coma, and Abell Clusters for our super- 
nova search, but in a few exceptional cases we chose to 
target specific galaxies (e.g. M31). In the present work, 
we remove this potential source of bias by selecting only 
the subset of SNe la found in the background of each 
field (z > Zf + 4000 km s -1 , where Zf is the redshift of 
the object targeted in a given field). During the search 
period, we detected 76 supernovae in total, and of these 
46 were classified as SNe la. After removing foreground 
and cluster supernovae, we are left with the 23 SNe la 
background events listed in Table [TJ Figures Q] and [D 
show our sky coverage with the locations of the discov- 
eries marked. Fields with the highest concentrations of 
nearby galaxies were typically observed every night as 
weather and season allowed. Over the first few years we 
observed the less rich fields on alternating nights. Ad- 
ditional time was allocated for the survey in the RSVP 
era, and fields were added at that time. 

For both TSS and RSVP, we employed image sub- 
traction techniques to reveal time vari able events. The 
TSS e mployed a modified version of the IPerlmutter et al.l 
(1999) sea rch code, and the RSVP used the process de- 
scribed in lYuan & Akerlofl (|2008f) . For our survey, we 
typically observe each field multiple times each night. 
From these we create 3 co-added frames representing the 
average of the first half of the nights data (NEW1), the 
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TABLE 1 

Background SN Ia Discovered with ROTSE-IIIb by Feb. 1, 2009 



IAU Name 


Disc. Date 


RA 


DEC 


z 


Disc. 
Mag 


^det 


Abs. 
Mag 


Used? 
LOSS 


Used? 

nrinn tt 

SDSd-11 


Note 


2004gu 


Dec 13, 2004 


12:46:24.7 


+ 11:56:56 


0.05 


17.4 


6 


— 19.4 


y 


y 






Mar 28 2005 


12:17:17.1 


+16:22:17 


0.02 


17.1 


10 


— 19.0 


y 


y 




2005ck 


Jun 05. 2005 


13:02:18.7 


+28:20:44 


0.09 


18.6 


2 


— 19.9 


n 


n 


a 


2005cr 


Jun 24^ 2005 


12:22:17.1 


+ 12:23:49 


0.02 


16.1 


10 


— 19.1 


y 


y 


c 


ZUUOllJ 


Ort 9005 

V7CL JU, ZUUO 




— 01T4T7 


u.uu 


1 7 fi 


I 7 


"1 7 
— iy. / 


y 


y 


b 


2005ir 


at r\r) / 1 / 1 / 1 — 

Nov 03, 2005 


01:16:43.7 


+00:47:40 


0.08 


18.5 


6 


— 19.4 


y 


y 


t 


2006an 


bco 21, zvOb 


1 rv -| A o o ""7 


+ 12:13:47 


0.06 


18.0 


2 


— 19.4 


n 


n 




2006cj 


May 17, 2006 


12:59:24.5 


+28:20:51 


0.07 


17.9 


11 


— 19.5 


y 


y 




2006ct 


May 25, 2006 


12:09:56.8 


+47:05:45 


0.03 


17.5 


11 


— 18.6 


y 


y 




2007if 


Aug 16, 2007 


01:10:51.4 


+15:27:40 


0.07 


19.5 


12 


-20.6 


y 


a 




2007kh 


Sep 07, 2007 


03:15:12.1 


+43:10:13 


0.05 


18.9 


10 


-19.3 


y 


y 




2007op 


TVT r\ A / 1 / 1 / 1 "7 

INov 04, 2007 


01:53:12.4 


+33:44:34 


0.09 


18.5 


6 


— 19.1 


y 


y 




2007qc 


Oct 27, 2007 


11:57:04.7 


+53:29:55 


0.04 


16.8 


11 


-19.1 


y 


y 




2007sp 


Nov 14, 2007 


12:04:42.3 


+49:11:09 


0.02 


16.8 


7 


-18.2 


y 


y 


e 


2007sw 


Dec 29, 2007 


12:13:36.9 


+46:29:36 


0.03 


16.0 


2 


-18.8 


ii 


n 




2008E 


Jan 04, 2008 


11:25:37.0 


+52:08:26 


0.03 


18.2 


18 


-18.7 


y 


y 




2008ab 


Jan 30, 2008 


11:34:45.9 


+53:57:51 


0.07 


18.2 


8 


-19.7 


y 


n 




2008ac 


Jan 30, 2008 


11:53:45.2 


+48:25:22 


0.05 


17.6 


15 


-19.2 


y 


y 




2008ad 


Jan 30, 2008 


12:49:37.2 


+28:19:47 


0.05 


18.0 


2 


-19.1 


ii 


n 




2008ar 


Feb 27, 2008 


12:24:37.9 


+ 10:50:17 


0.03 


16.9 


17 


-18.8 


y 


y 




2008bg 


Mar 12, 2008 


12:51:11.9 


+26:17:40 


0.06 


18.7 


9 


-19.3 


y 


y 


c 


2008by 


Apr 19, 2008 


12:05:21.0 


+40:56:46 


0.05 


17.2 


2 


-19.8 


ii 


n 




2008bz 


Apr 22, 2008 


12:38:57.7 


+11:07:46 


0.06 


17.7 


6 


-19.2 


y 


y 





NOTE. — "Also reported by LOSS (JPugh ct al. 2005). 6 Also detected by SDSS -II and used in their volumetric rate 
calculations ( Dilday cl al. 2008, 2010). c Spcctroscopically confirmed by the CfA (Modiaz ct al. 2005: Colcsanti ct al. 
120061 : lYuan et al.|[2008bl) . d Also reported by the Nearby SN Factory I IYuan et al JI2007bt TScalzo et alJI2010l) . 'Maximum 
light not well constrained. Magnitudes given arc in the ROTSE-III bandpass calibrated against the USNO-B1.0 R2 
mags as described in the text. The AT det column gives the number of times each object was selected by the automated 
search pipeline. The ninth and tenth columns indicate whether the object was included in the rate calculations using 
the LOSS pAMD or the SDSS-II model pAMD, respectively. 

limiting magnitude for the nightly coadds was typically 
around 18.3 mag. 

We further vet potential SN candidates against cata- 
logs of known time- variable phenomengQ. Most of our 
fields were observed by the SDSS as well, so we were 
able to remove contamination from un-cataloged vari- 
able stars by checking for matching "host" stars in the 
SDSS data, which is complete to ~ 4 magnitudes be- 
low our survey depth. Remaining sources were vetted 
by human scanners. All sources judged as possible tran- 
sients were spectroscopically classified by us or others in 
the community. This is a key distinguishing feature of 
our survey; other wide-field transient searches generate 
more candidate supernovae than can be classified with 
the spectroscopic resources available, so a filtering pro- 
cess such as light curve fitting is typically employed to 
select only a subset for final spectroscopic confirmation. 

Table Q] also notes the number of times that each dis- 
covery was recovered by the automated search pipeline 
(Ndct)- A final analysis of the data adds low signifi- 
cance (i.e. < 5cr) detections that are not passed by the 
automated selection cuts, and we determine the peak 
magnitudes for each supernova in our sample by fitting 
light curve templates to these final data sets. We use 
R-band light curve templates for SNe Ia from LOSS, 
which were constructed by interpolating between the 
faint/fast SN 1 999by and the bright/slow SN 1991T 
(|Li et alllmlbl ). See g3] for a discussion of our unfil- 
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Fig. 2. — Similar to Figure [T] but for the Summer- Fall fields. 

average of the second half (NEW2), and the average for 
the whole night (NEW). We use filtering to remove par- 
ticle events and other artifacts during the co-addition. 
We then subtract the properly convolved reference tem- 
plate from each of these to generate subtracted frames 
(SUB1, SUB2, and SUB). To remove contamination from 
solar system bodies and remaining artifacts, we require a 
candidate to be detected (^ 2.5a) at consistent celestial 
coordinates on each of SUB1 and SUB2. The candidate is 
finally required to be detected at the 5cr confidence level 
on the combined subtraction, SUB. The typical FWHM 
for the survey was 3.2 pixels (about 10") and the 3a 



7 For example: http://scully.cfa.harvard.edu/cgi-bin/checkmp.cgi 
http: / /simbad.u-strasbg . f r/simbad/sim-f id 
http : //nedwww . ipac . caltech.edu/ 
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tered spectral response and a comparison to the R-band. 
We use mpf it . pro in IDL to determine the best fitting 
template shap e, and subtract off th e Galactic extinction 
in the R-band f Schlegc l et al.lll998l ). We employ the dis- 
tance modulus (based on independent distance indicators 
when available through NEE0 or otherwise based on the 
redshift in a flat, H = 71 kms -1 Mpc -1 , Q m = 0.27 
cosmology) to determine the peak absolute magnitudes 
shown in table [TJ SN 2004gu was first observed near 
peak, which adds some uncertainty to the date and mag- 
nitude of maximum light; however, spectroscopic phase 
information is consistent with the derived date of maxi- 
mum light, and thus we expect the peak magnitude esti- 
mate to be accurate. According to its spectroscopic age, 
another event, SN2007sp, was first detected about two 
months after maximum light, so the light curve fit likely 
underestimates the peak. 

2.1. Sample Notes 

Although our sample is small, it contains several SNe la 
worth of special note. 

• SN 2004gu is a significant outlier on the Hub- 
ble diagram, and other authors have grouped it 
with peculiar events like SN 2006gz and SN 1999ac, 
which may or may not be linked to super- 
Chandrashekar explo sions (jContreras et al.l 120101 : 
iSilverman et alJ[20lTh . 

• SN 2006ct showed the characteristic features 
of a Type la supernova, but the blue-shift of 
the line minim a near maximum lig ht was only 
- 6000 km s^ 1 (|Quimbv et al.l I2006T) . The spec- 
tra show s ome similarities to the hig hly peculiar 
SN 2002es (jGaneshalingam et al.ll2012l ). which was 
also hosted by an early type gal axy, but the p ho- 
tometric decline of SN 2006ct (|Quimbvl I2006T) is 
slowe r and more remini scent of SN 2002cx (|Li et al] 
l200l and SN 2005hk (IPhillips et al.ll2007l) . 

• SN 2007if is the most luminous Type la super- 
nova known, and it may be the result of a super- 
Chandrashekar m ass explosion (Sc alzo et al.ll2010l : 
lYuan et alJl201Cft . 

• SN 2007qc was hosted by an extremely low lumi- 
nosity dwarf galaxy (see At Mb ~ — H, the 
host is among the faintest detected for any SNe la. 

• SN 2008ar showed spectroscopic evidence for high 
velocity material and u nusually strong ab sorption 
in the Ca II IR triplet (lYuan et al.l I2008al^ . w hich 
is reminiscent of SN 20071e (jSimon et al.ll2009h . 

In addition, at least 7 of our SNe la (2004gu, 2005ck, 
2005hj, 2006cj, 2007if, 2008ab, 2008by and possibly 
2005bg) show cither unusually high luminosities or are 
spectroscopically classified as SN 1991T/1999aa-like. 
This appears to be a high fraction even for a magni- 
tude limited survey when com pared to the LO SS census 
of SNe la in targeted galaxies (|Li et al.ll2011b|h 

8 http://ned.ipac.caltech.edu/ 



3. SN la PSEUDO-ABSOLUTE MAGNITUDE 
DISTRIBUTIONS 

Our survey is flux limited, and as a result we have a 
bias against selecting low luminosity SNe la. This bias 
is exacerbated by our removal of discoveries made at or 
below the redshift of the targeted galaxy clusters. To 
determine the event rate for the SNe la population as a 
whole, we must therefore correct the number of events in 
our observed sample for the fraction we are likely to have 
missed due to, among other concerns, their lower lumi- 
nosities (see fj5]) . This requires knowledge of the intrinsic 
luminosity function for all SNe la in our search volume 
combined with the host galaxy absorption distribution. 

In this section, we consider the pseudo- absolute magni- 
tude distribution (pAMD) for SNe la, which is the distri- 
bution one would obtain by correcting the peak observed 
magnitudes of a complete sample of SNe la for distance 
and Galactic extinction but not for host galaxy extinc- 
tion. This distribution can directly be used to test the 
selection efficiency of a flux limited survey (see §5§. First, 
we consider in §3.11 how the available SNe la luminosity 
functions derived (mostly) through filtered observations 
may compare to our unfiltered survey data. We then dis- 
cuss the pAMD compiled b y LOSS from th eir targe ted, 
but volume limited sample (|Li et al.ll2011bfj in §TL2j and 
we derive an empirical pAMD from a volume limited sub 
sample of the SDSS-II survey in §3.31 We find that the 
LOSS and SDSS-II pAMDs are not consistent, which is 
discussed in H3A[ 

3.1. Comparison of R-band and ROTSE-IIIb Unfiltered 
Magnitudes 

We calibrate our instrumental magnitudes against the 
USNO-B1.0 R2 system since some of our search fields 
are outside of the SDSS footprint. For fields covered 
by both the SDSS and USNO-B1.0, we compared the 
USNO-B1.0 photometry to SDSS photometry converted^ 
to the R-band using R = r - 0.2936(r - i) - 0.1439, 
and we find the two systems agree to about 0.01 mag 
on average and have a field to field dispersion of about 
0.14 mag. Thus our unfiltered ROTSE-IIIb magnitudes 
should be, on average, similar to the R-band, but objects 
with spectral energy distributions that differ significantly 
from the field stars employed in the calibration may be 
offset from their true R-band values. 

We estimate the potential offset between our unfiltered 
magnitudes and the true R-band system for SNe la as fol- 
lows. We adopt the spectral energy distribution of a G5V 
star (|Picklesl 1T998) as representative of a typical field 
standard used in calibrating the photometric zeropoint 
of our unfiltered imaging data. We scale this template to 
zero magnitude in the R-band, and then measure the syn- 
thetic flux using the approximate tran smission function 
of ROTSE-IIIb (jQuimbv et al.ll2007bn . which is plotted 
in Figure 03 We next scale a template^ spectrum for a 
Branch-normal SN la at maximum light to zero magni- 
tude in the R-band and measure its flux in the ROTSE- 
IIIb bandpass. The ratio of this and the scaled field star 
flux (also measured in the ROTSE-IIIb bandpass) gives 
the expected departure of the ROTSE-IIIb magnitudes 

9 See the Lupton (2005) equations at 

http: / /www. sdss . org/dr7/algorithms/ sdssUBVRITransf orm .html 
lu http : // supernova, lbl .gov/~ nugent/nugent_templates .html 
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3000 4000 5000 6000 7000 8000 9000 10000 
Wavelength (A) 

Fig. 3. — Determination of the magnitude offset between filtered 
and unfiltered measurements. The top panel shows the spectral 
energy distributions of a typical field standard (blue) compared 
to a high luminosity SN la (purple) and a low-luminosity SN la 
(orange). These can be convolved with the ROTSE-IIIb (grey) 
and R-band (dashed red) response curves in the lower panel to 
determine the magnitude offset. The dotted green curve is the 
typical response for an Apogee7 camera, which is used by KAIT 
for the LOSS. See text for details. 

from the true R-band system. We measure this offset 
in the range 0.02 < z < 0.09 by redshifting the SN la 
template, and find that our unfiltered magnitudes may 
be 0.05 to 0.10 mag brighter than the true R-band sys- 
tem. We rep eat this process for a high luminosity SN la 
(SN 1999aa; lGaravini etlrtl 120041), and a low luminosity 
event (SN 1991bg; lFilippenko et aLlll992t) . We find that 
a SN la with the same R-band brightness as a field stan- 
dard will generate a magnitude in the unfiltered system 
(calibrated against the R-band) that is within ^0.1 mag 
of its true R-band magnitude. 

While this shows that our observed magnitudes should 
be similar to the true R-band system, a more pertinent 
concern is how the derived absolute magnitudes will com- 
pare. In particular, we have calculated absolute magni- 
tudes from our unfiltered data by simply subtracting off 
the distance moduli and Galactic (R-band) extinction 
terms; we have not included a "k-correction" term to ac- 
count for the changing rest frame bandpass of our unfil- 
tered system over the modest rcdshift distribution of our 
sample. Below, we compare the pseudo absolute magni- 
tudes derived from our data to rest frame (k-corrected) 
R-band magnitudes calculated for the same sources based 
on filtered photometry to check for an o ffset. 

The CfA3 sample includes MLCS2k2 (iJha et al.ll2007l) 
fits fo r four of the supernovae in our sample (jHicken et al.l 
I2009D. including two that were also in the SDSS-II sample 
(jKessler et al.1 [20091 ) . We calculate the pseudo- absolute 
magnitudes from these fits using the A parameters and 
adding back in the host extinction term. For this compar- 
ison, we calibrate the ROTSE-IIIb magnitudes against 
the SDSS field stars (converted to the R-band as de- 
scribed above) to remove the USNO-B1.0 vs. SDSS-II 
offset, which we find to be biased in this small overlap- 
ping sample. 

Our pseudo absolute magnitude estimate for SN 2005hj 
is about 0.1 mag brighter than the fits from the CfA3 and 
the SDSS-II. For SN 2005ir, the SDSS-II and CfA3 differ 



from each other by about 0.1 mag, which is compara- 
ble to the combined measurement error. Our absolute 
magnitude for SN 2005ir is about 0.02 mag brighter than 
the SDSS-II value and 0.09 mag fainter than the CfA3 
value. Our measurements of SN 2006an and SN 2006cj 
are 0.03 mag brighter and 0.16 mag fainter than the CfA3 
peak magnitudes, respectively. 

Additionally, the Nearby Supernova Factory has re- 
ported a peak V-band magnitude and V — R color for 
SN 2007if that implies a peak absolute R-band maxi- 
mum of about —20.4 mag. Calibrating again st the SDSS 
r-band AB magnitudes, lYuan et al.l (|2010T ) found the 
ROTSE-III observations peaked at M r = -20.4. Field 
stars will typically "brighten" by about 0.22 mag when 
converting from the SDSS r-band system to Cousins R- 
band (Vega) magnitudes, which shifts the target pho- 
tometry by the same amount. Our estimate for the peak 
of SN 2007if is approximately 0.15 mag brighter than the 
Nearby Supernova Factory value. 

To, summarize, of the 7 pseudo-absolute magnitude 
measurements available in the literature based on (k- 
corrected) filtered photometry, our ROTSE-IIIb esti- 
mates are ~ 0.1 mag brighter in three cases, ~ 0.1 mag 
fainter in 2 cases, and roughly the same in the remaining 
2 cases. On average, our values are 0.02 mag brighter. 
We conclude that there is no significant systematic off- 
set between our unfiltered, pseudo-absolute magnitude 
estimates and the (mostly) R-band system that defines 
the luminosity functions considered below ( ^3.2l and l3.3|l . 
and these should agree to about ~ 0.1 mag for individ- 
ual events. We address the impact the possible system- 
atic error in our magnitude system with respect to the 
LOSS and SDSS-II pseudo-absolute magnitude distribu- 
tions will have on our rates in [JH 

3.2. The LOSS pseudo-Absolute Magnitude Distribution 

for SNe la 

The Lick Observatory Supernova Search (LOSS) has 
recently published a distribution of pseudo-absolute mag- 
nitudes for SNe la from their volume limited search 
(|Li et al.l f2TTTTbh . Like ROTSE-IIIb, the LOSS survey 
engine (the Katzman Automatic Imaging Telescope, or 
KAIT) runs with an unfiltered CCD imager, but follow- 
up observations are often obtained through standard 
broadban d filters. The S Ne la luminosity function re- 
ported bv lLi et all (|2011bl ) is primarily constructed from 
R-band data, which the authors describe as the best 
match to their unfiltered survey data, and also from unfil- 
tered observations as well. In the lower panel of Figure [3] 
we show the response curve for the Cousins R-band com- 
pared to the typical unfiltered response of an Apogee7 
camera, which is used by LOSS. 

It is important to note that the LOSS SNe la luminos- 
ity function was drawn from a sample that was highly bi- 
ased with respect to the host galaxy properties. In partic- 
ular, the LOSS targeted mainly high luminosity galaxies. 
It has been demonstrated t hat certain SNe la su b-types 
prefer certain host galaxies (jSullivan et al.l 120061 ). so the 
luminosity distribution of the LOSS sample may be bi- 
ased with respect to the larger population from which 
our sample is drawn. 

An obvious difference between the LOSS and ROTSE- 
IIIb SN la samples is the presence of high luminosity 
{Mr < —19.6) events in the latter. Although we could 
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follow previous works and simply discard our most lumi- 
nous SNe la, we instead chose to calculate an inclusive 
rate. We must, therefore, augment the LOSS luminosity 
function to account for such events or else our effective 
survey volume will be under estimated. Lacking detailed 
demographics for this population, we simply assume that 
the population of high luminosity SNe la represents only 
about 1% of SNe la in a given volume; if it were much 
more these would be detected more frequently, and if it 
were much less we could not expect to find any in our 
small sample (see $5|. We assume that SN 2007if was 
a particularly luminous example of this population (as 
noted above, it is the most luminous SN la known) and 
choose a half Gaussian with a peak at Mr = —19.6 and 
a = 0.4 mag to stand in for the unknown luminosity dis- 
tribution. We will use this augmented version of the 
LOSS pAMD to calculate the SNe la rate in S0 

3.3. The SDSS-II pseudo-Absolute Magnitude 
Distribution for SNe la 

For our rate calculations in fJTl we will also use the 
same bimodal-Gaussian luminosity function attenuated 
by a host galaxy absorption distribution that was as- 
sumed by IDildav et all (|2010D in deriving the SDSS-II 
rate. The mode l is a revised yersiorO of the distribu- 
tion assumed by IDildav et all (120081) th at takes into ac- 
count the findings of lKessler et all (|2009() . The difference 
makes for only a negligible change in our rate measur- 
ment. The intrinsic SNe la luminosity function is defined 
in terms of the MLCS2k2 light-curve shape/luminosity 
parameter, A. The distribution of A is defined as a 
Gaussian with a = 0.19 for A < -0.2 and a = 0.40 
for A > —0.2, and it is truncated to lie within the range 
(—0.4 < A < 1.8). The final, pseudo- absolute magnitude 
distribution is achieved by adding host galaxy absorp- 
tion drawn from a distribution with P(Ar) oc q-Ar/oss^ 
With our choice of Ho , peak pseudo-abs olute magnitudes 
are then defined as (cf. Uha et aLll2007l ): 

Mr = -19.313 + 0.579A + 0.254A 2 + A R (1) 

When applying the SDSS-II pAMD model to our 
ROTSE-IIIb discoveries, it is important to note that 
the SDSS-II rate is based on a sub-sample of SNe la, 
and to be consistent with the SDSS-II pAMD model 
we must also remove objects from our sample. In par- 
ticular, the SDSS-II includes only the fraction of the 
population that is well fit by MLCS2k2 (V &t > 0.001) 
with a light curve width parameter A > —0.4. Such 
a cut excludes high lumin osity SNe la, and i n particu- 
lar, SN 1999aa-like events (jKessler et alll2009f ). We note 
that this is a common choice. The Supernova Legacy 
Survey (SNLS), for example, has presented a volumetric 
rate for SNe la at z ~ 0.45, but the spectroscopically 
confirm ed Type la SN 2003fg (SNLS-03D3 bb) is not in- 
cluded (jNeill et al.ll2006l : Howell et alll2006l ). Indeed, the 
selection cuts that are applied in selecting candidates for 
spectroscopic follow-up may also pose a bias against high 
luminosity SNe la. We therefore remove events with peak 
luminosities brighter that Mr < —19.7 when using the 
model pseudo-absolute magnitude distribution from the 
SDSS-II. 



We next look at the actual distribution of pseudo- 
absolute magnitudes from a volume limite d sub-sample of 
the first year SDSS-II SNe la discoveries (Fricma n et all 
120081: ISako et~al1l20r)8h . which can be compared to the 
as sumed model as we ll as the LOSS pAMD. According 
to IDildav et"all (|2008| ). the SDSS-II could have detected 
even a sub-luminous (SN 1991bg-like) SN la out to a red- 
shift of z ~ 0.12, so we assume the survey to be volume 
limited out to this limit. There were 29 spectroscopically 
confirmed (or probable) SNe la from the first year of the 
SDSS -II sear c h wit hin this redshift limit. In addition, 
IDildav et all (l2008h lists one probable SN la (SDSS-II 
SN09266; based on photometric constraints) with a host 
redshift of z = 0.0361. This so urce is not included in the 
sample of Dil dav et all (|2010T ) due to an additional cut 
on the photometric screening. A second photometrically 
probable source, SDSS-II SN09739, with a spectroscopic 
redshift for its host and a third, photometrically proba- 
ble SN la (SDSS-II SN11092) are also given for a total 
of 32 SNe la in the volume limited sample. 

We now derive an empirical SNe la pAMD from this 
volume limited sub sample of the first-y ear SDSS-II. 
We ad opt the temporal s electi on cuts of Dilda v~et all 
(2008) and iKessler et al.l ([2009) to remove events with 
poorly constrained maxima, but unlike these samples, 
we include peculiar events. For the bulk of the sam- 
ple, we derive the pe ak absolu t e R-b and magnitudes 
from the MLCS2k2 (Uha et all 12007ft fit parameters 
given in IKessler et al.l (|2009l ). We add the host ex- 
tinction estimate back in to retain the pseudo-absolute 
distribution required for rate studies. Two peculiar 
SNeJa J _SN 2005hk and SN 2005gj, were cut from the 
IKessler et all (|2009f) sample; it is still a matter of de- 
bate if these events truly belong to the Type la class 
(e.g iValenti et al.l 120091: iFqlev et all 12010 1 : iMaund et all 
1201ft lAldering et all 120061 : iTrundle et all 120081 ). We in- 
clude them here for completeness. For SN 2005hk, we 
take the L OSS measureme nt for the R-band absolute 
magnitude (|Li et alll20TTbT ). and for SN 20 05g.j we take 
the pe ak absolute r-band magnitude from iPrieto et all 
(2007) and convert this to R-band. We similarly use the 
SDSS- II photometric measurements of iHoltzma n et all 
( 2008) to estimate peak absolute R-band magnitudes for 
two spectroscopically probable SNe la: SN 2005je and 
SDSS-II SN06968. Finally, we choose to include an esti- 
mated peak magnitude for SDSS-II SN09266. There is no 
published light curve or peak magnitude for this event, 
so we assign it an ad hoc pseudo-absolute magnitude of 
Mr = —1 6 mag based on th e heavy extinction reported 
(A„ ~ 4: IDildav et al.ll2008h . 

We thus have peak magnitudes for 22 of the first 
year SDSS-II SNe la; eight of the remaining events 
have poorly constrained maxima. The portions of the 
light curves sampled by the SDSS-II, however, imply 
a similar distribution of peak magnitudes to the well 
constrained sample. Photometry is not available for 
the other two photometrically selected events (SDSS-II 
SN09739 and SN11092), so we cannot include these. The 
peak R-band pseudo-absolute magnitudes (corrected to 
H a = 71kms- 1 Mpc- 1 ) for the volume limited SDSS-II 
sample are listed in table [2 
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TABLE 2 

SDSS-II First Year Volume Limited 
Sample 



SDSS ID IAU Name Mr Note 



722 


900 ^nH 




739 


200 5of 




774 


2005ei 




1241 


2005ff 


-18.97 


1371 


2005m 


-19.39 


2102 


2005m 




2561 


2005fv 


-19.03 


3256 


2005hn 


-19.11 


3592 


2005gb 


-19.30 


3901 


2005ho 


-19.33 


4524 


2005gj 


-20.47 


5395 


2005hr 


-19.36 


5549 


2005hx 


-19.23 


5944 


2005hc 


-19.41 


6057 


2005if 


-18.95 


6295 


2005js 


-17.72 


6558 


2005hj 


-19.36 


6773 


2005m 


-19.28 


6962 


2005jc 


-19.28 


6968 




-19.06 


7147 


2005jh 


-18.96 


7876 


2005ir 


-19.29 


8151 


2005hk 


-18.36 


8719 


2005kp 


-19.30 


9266 




-16.00 


9739 






10028 


2005kt 


-19.06 


10096 


20051j 




10434 


20051k 




10805 






11067 


2005ml 




11092 







Note. — a Peak magni tude con- 
verted from iPrieto et all J2007D. b Peak 
magnitude estimated from Holtzman ct al. 
(2008) light curve. c Pcak magnitude from 
I Li et al.| j2011b|) . d Peak ma gnitude esti- 
mate based on description in Dilday ct al. 
j2008f). Remaini ng magnitude estimates 
from lKessler et al.l ( 120091') . Missing Mr val- 
ues result from events discovered at the be- 
ginning or end of the search period for which 
maximum light is not well constrained, or 
because photometry is not available (9739 
and 11092). Also note that the values in the 
Mr column arc pseudo-absolute magnitudes 
that have not been corrected for host galaxy 
extinction. 

3.4. Comparison of the LOSS and SDSS-II 
pseudo- Absolute Magnitude Distributions 

A comparison of the LOSS SNe la luminosity function 
and the SDSS-II first year volume limited sample shows a 
striking result: the samples do not agree (Fig. [4j) . About 
66% of the LOSS SNe la are fainter than —19 mag, but 
only about 27% of the events in the SDSS-II sample are 
this faint. We discuss the implications of this disagree- 
ment in Only SN 2005gj, which appears to draw 
its power through interaction with a hydrogen shell and 
could technically be classified as a Typelln, approaches 
the lumino sity of our higher luminosity events , such as 
SN 2007if (jScalzo et al]l20Tot lYuan et al]l20l0h . 

Also shown in Figure |4] is the R-band pAMD model 
that was assumed in calculating the SDSS-II SNe la rates 
HDildav et all 120101) . which also appears to differ from 
the selected SDSS-II population. From a KS test, we 
find that there is only about a 4% chance that the vol- 




-16 -17 -18 -19 -20 -21 

SN la Peak Pseudo-Absolute Magnitude 



Fig. 4. — Absolute SNe la magnitude distributions from the 
LOSS (blue) and SDSS-II (red) volume limited samples. The pink 
line is the subset from the volum e limited SDSS-II s ample that 
was used in the rate calculation of Dilday ct al. (2010). The dot- 
ted orange line is the model pAMD assumed in the SDSS-II rate 
calculation; we use this model to determine a SNe la rate that 
can be directly compared to the SDSS-II value. For the pseudo- 
absolute magnitudes, the observed magnitudes are corrected for 
Galactic extinction and the distance moduli are removed, but the 
host extinction is left uncorrected. 

ume limited SDSS-II SNe la sample is co nsistent with 
the at tenuated, bimodal-Gaussian model of iDildav et al.l 
(2010). The reddening distribution assumed in the model 
predicts that it is extremely improbable to have highly 
extincted events such as SDSS-II SN09266 (again, this 
was not spectrosco pically confirmed and not included in 
IDildav et al.ll2010T) . We repeated the KS test using only 
the member s of the volume lim ited sample that were 
also used by IDildav et al.l (|2010D in their rate calculation 
and find the probability that these events were drawn 
from the assumed model is only about 3%. Except for a 
few low-luminosity events, the observed SDSS-II sample 
is systematically biased to higher luminosities than the 
model predicts. We discuss the SDSS-II model as applied 
to the ROTSE-IIIb sample in gH 

4. DETECTION EFFICIENCY 

In this section we determine the detectability of sources 
as a function of observed magnitude relative to the lim- 
iting magnitude of a given image. This curve will be 
used later to determine the probability of detecting sim- 
ulated sources with various observed magnitudes, which 
is needed to calculate our overall survey efficiency. 

To determine the probability that a source of a given 
brightness will be detected, we perform a Monte Carlo 
simulation. We add simulated point sources to our co- 
added images and determine the fractions recovered as 
a function of magnitude. This is done relative to the 
limiting magnitude of the images and coadds recorded 
in the logs, which is calculated from the 90th per- 
cent ile magnitude of all obj ects extracted by SExtrac- 
tor (|Bertin fc Arnouts1ll996l ). Our reference images are 
typically constructed from dozens of the best individ- 
ual frames from each field, so they contribute negligi- 
bly to the noise on the subtracted frames. We therefore 
choose to calculate the detectability of sources directly 
from blank sky regions of the co-added images without 
performing the computationally expensive process of im- 
age subtraction. A consequence of this is that we do not 
account for the effects of contaminating host galaxy light. 
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Other surveys that employ image subtraction have previ- 
ously tested the effects of host light on detectability and 
found very little change from bright hosts to blank sky 
(jNeill et al.ll2006[ ). Nonetheless, our procedure may over- 
estimate our detection probability for supernovae hosted 
by bright galaxies (i.e. the rate may be underestimated). 
Our results do, however, have direct applications for ob- 
jects with faint hosts or otherwise dark backgrounds. 

On a given image, the sensitivity to objects of a given 
magnitude will vary with location in the field due to vari- 
ations in the instrumental PSF as well as structure in the 
atmospheric transmissivity (e.g. passing clouds), which 
can be resolved by our wide field. To account for these 
effects, we locate a set of isolated field standards on each 
image to use in calculating the local PSF and zeropoint. 
We divide each image into a grid with as many cells as 
can be made while maintaining around 100 or more field 
standards in each (typically 3x3 cells), and then cal- 
culate the average zeropoint and PSF for each of these 
cells. The PSF is found using the DAOPHOT routines in 
IDL, and as noted above f ffTTj) . we use the USNO-B1.0 
R2 system for our zeropoints. 

To locate appropriate blank sky regions at which we 
may add our simulated point sources, we first run SEx- 
tractor on each unadulterated image with a low detection 
threshold to identify all objects and artifacts in the field. 
We then convolve the object mask image from SExtractor 
with the average PSF and select from the remaining, un- 
masked pixels when placing simulated sources. We find 
that even the small perturbations induced by adding a 
very faint (m > my lm + 2) test source can sometimes el- 
evate a faint (m ~ mij m ) source previously ignored by 
SExtractor into a spurious detection. It is therefore im- 
portant to avoid locating test sources near such areas, 
which is why we use a lower extraction threshold when 
creating the object mask. 

We use the local PSF and zeropoint determined above 
to determine the shape of the test source and scale it 
to a randomly chosen magnitude. Using the dao_value 
routine, we create a thumbnail image of the source. We 
then add this source to the image one pixel at a time 
by drawing from a Poissonian distribution with the ex- 
pectation values set by the thumbnail pixel values and 
accounting for the detector gain. We finally run SEx- 
tractor with the same parameter file used in the search. 
In order to simulate a statistically meaningful number of 
test objects without biasing the extraction process, we 
place only 20 test sources in each cell and repeat the 
simulation multiple times for each image. The minimum 
pitch for the test objects is set to be at least four times 
the average FWHM, so they do not interfere with each 
other. For the survey, candidates must be detected at a 
signal to noise greater than a set limit and they must have 
a FWHM between half and double the image average to 
be considered. Both of these quantities are effected by 
the noise in the image and Poisson statistics. We there- 
fore define our detection efficiency from the number of 
simulated point sources that are detected by SExtractor 
with parameters that pass these selection cuts relative to 
the total number simulated. We performed this test on 
over 7000 co-added images made from 2 or more frames. 
For comparison, our survey includes about 55000 master 
subtractions. 

Figure [5] shows the average 5er detection efficiency 




Magnitude Relative to m. 



Fig. 5. — Detection efficiency for simulated point sources as a 
function of magnitude relative to the logged limiting magnitude, 
m lim- The top panel shows that the shape of the detection effi- 
ciency curve is nearly independent of the limiting magnitude for 
a fixed FWHM. The lower panel shows that for a fixed mi; m , the 
shape of the detection efficiency curve depends strongly on the av- 
erage FWHM of the image. For both plots, the 50% completeness 
is brighter than the logged limit mainly because the simulations 
require a 5(7 detection while the logged limits are derived from all 
extracted sources (i.e. 3<r detections are included). 

curves for isolated objects. As can be seen, the shape 
of the curve is robust to changes in the limiting mag- 
nitude, but the effect from variations in the FWHM is 
quite noticeable. 

Accurate limiting magnitudes were not recorded for the 
co-added images during the first half of our survey, but 
we do have the limits for the individual images. Using 
our later data, we determine an empirical relation be- 
tween the limits of the individual images and the coadds 
derived from them. We started by assuming the limiting 
magnitudes for each individual image is approximately 
equal to three times the sky noise in the aperture. In- 
verting this, we can find the sky noise expected on the 
co-added frame from the sum of the individual images, 
and turn this into a limiting magnitude guess. Compar- 
ing these guesses to the actual values recorded in the 
logs over the second half of our survey, we find that the 
agreement is good for the deepest data, but there is a 
systematic offset for more shallow coadds. We perform 
a linear fit to these data and use it to derive magnitude 
limits for the coadds created in the beginning of the sur- 
vey. 

5. SURVEY EFFICIENCY 

It is not enough to merely detect a source to count it 
as a discovery; each such candidate must pass additional 
screening to be included in our sample. This includes 
both machine cuts designed to reject false positives and 
cuts made by the humans who have the final say in elevat- 
ing candidates to discoveries. Often it is straight forward 
to determine the effects of the machine cuts through a 
Monte Carlo process: simulated supernovae drawn from 
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the expected population are "observed" using the actual 
survey cadence and limiting magnitudes, and the result- 
ing "candidates" are passed through the same machine 
vetting process as the real data. The sensitivity of the 
survey to sources with various apparent brightnesses and 
light curves can then be directly determined by noting 
the fraction of the simulated supernova that pass. The 
effects humans play in the ultimate selection of candi- 
dates could in principle be determined through a similar 
simulation if the test were performed concurrently with 
the actual survey to account for learned behavior. To 
get statistically meaningful results, the human scanner 
would need to be presented with an enormous volume 
of simulated data, which is far too large of a burden on 
the researchers to be done in practice, although previous 
studies have performed limited versions of such human 
experimentation, which guide our tests. 

As stated above, the basic machine cuts applied to the 
raw candidates are 1) the FWHM should be between 0.5 
and 2 times the image average, and 2) the signal-to-noise 
should be at least 2.5a on SUB1 and SUB2, and 5er for 
SUB. To help reduce the number of false positives created 
by imperfect subtraction of galaxy light, we also required 
that candidates near the cores of cataloged galaxies show 
an increase of at least 5% in flux over the reference tem- 
plate (given the bright limits of our survey and prior 
coverage of our search fields by galaxy surveys, it is fair 
to assume that almost all galaxies as bright or brighter 
than our SNe la sample were cataloged). We are mainly 
sensitive to SNe la brighter than —19.0 absolute, so hosts 
fainter than about —20.3 are unaffected by this cut. Of 
course, when the host galaxy is resolved, only the frac- 
tion of light coincident with the candidate matters in this 
percentage increase cut. For example, we recorded a 50% 
increase for SN 2005bg (M pca k = —19), which was dis- 
covered near the core of a M r ~ —21 host. Candidates 
coincident with sources not listed in the galaxy catalogs 
(typically foreground stars) were required to show a 15% 
increase. The remaining candidates located more toward 
the outskirts of cataloged galaxies or in isolated regions 
were not required to show a minimum increase in flux. 
We do not include this selection criteria in our Monte 
Carlo simulations, so our efficiency in discovering super- 
novae in luminous host galaxies may be overestimated, 
which again means that our rate in giants may be under- 
estimated. 

The basic setup for our efficiency calcul ation is simi- 
lar to that used in previous rate studies (cf. IDildav et al.l 
20j08|. We determine the efficiency as a function of lu- 
minosity distance by considering a series of thin shells 
over which the efficiency (and rate) can be assumed to 
be constant. In each of these shells we simulate a num- 
ber of SNe la randomly drawn from the pseudo- absolute 
magnitude distributions discussed in $3l choose an ap- 
propriate light curve template and scale it to the given 
distance, and select a random date for maximum light. 
By comparing the expected magnitudes on the dates 
each simulated supernova would have been observed by 
our survey and using the actual survey detection limits 
on those nights, we determine the fraction of simulated 
SNe la that could have been detected in each distance 
bin. We take into account the detection efficiency curve 
discussed in Sj4] appropriate to a given observation based 
on the limiting magnitude and FWHM recorded (or es- 



timated). For each observation epoch of each simulated 
supernova, we draw random numbers and compare these 
to the detection probabilities to determine if a simulated 
detection is made on each of the SUB1, SUB2, and SUB. 
For each of our 177 fields, we simulated 100 SNe la in 
each of 360 distance shells spaced logarithmically from 
40 to 1000 Mpc, which totals to over six million simula- 
tions for each pAMD assumed. To naturally account for 
the changing number of survey fields, we allow simulated 
supernovae to peak anytime in any field during the entire 
survey period including off season. 

To assess how accurately these simulations model the 
actual survey, we compare the luminosity distance, dis- 
covery magnitude, peak absolute magnitude, and detec- 
tion number distributions predicted by the Monte Carlo 
simulations to the observed distributions. Figure |6] shows 
that the survey model discussed thus far is a poor rep- 
resentation of the true sample if the augmented LOSS 
pAMD is adopted (the same result holds if the SDSS-II 
model pAMD is used). The model predicts the SNe la 
should be further away than observed, discovered at mag- 
nitudes fainter than observed, and the simulated super- 
novae are typically detected on fewer nights. This last re- 
sult is key to understanding the failure of the model thus 
far. The model predicts that > 35% of the SNe la should 
be detected by the automatic search selection cuts on 
exactly one night. After all, the largest volume element 
is for the most distant events that are only detectable 
at their peak. Since we have zero SNe la detected on 
just one night, it would seem that our survey is biased 
against such "one-nighters." A possible explanation for 
this is the influence of the human scanners who have the 
final say in selecting targets. Although it was never a for- 
mal requirement of the search that events be detected on 
more than one night (and in fact, we triggered spectro- 
scopic follow-up more than once on the same night based 
on a first detection), the informal requirement that can- 
didates "look good" may ha ve created a b i as ag ainst the 
weakest detections. Indeed, IDildav et al.l (|2008T ) studied 
the response of human scanners to simulated supernovae 
and determined that there is such a bias against one- 
nighters. 

To account for our apparent selection bias against weak 
detections, we apply a cut to the simulations and the real 
sample alike on the number of nights detected. We find 
that setting a minimum number of 5 nights, the model 
and actual distributions agree for both the augmented 
LOSS pAMD and the SDSS-II model. This cut prefer- 
entially removes the most distant supernovae with the 
faintest observed magnitudes, so after applying this sin- 
gle correction, all of the model distributions considered 
agree with the actual observed sample (Fig. [TJ) . The min- 
imum number of detections was set by calculating the 
probability that the observed sample was drawn from the 
model distribution via a KS test. We assign a consistency 
probability, Pks, by randomly drawing from the model 
distribution and comparing the maximum displacement 
in the cumulative distributions of the model and each 
draw. For 4 or fewer detections and the LOSS pAMD, 
only 0.3% (or 0.4% with the SDSS-II pAMD model) or 
less of the random draws have a displacement as large or 
larger than the real data compared to the models, so we 
reject these models. With a minimum of 5 detections, 
the probability rises to 9% (5%), so we cannot reject the 
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Fig. 6. — Distributions predicted from the Monte Carlo simulations (grey) compar ed t o the observ ed sample (gree n) before correction 
of selection bias. The augmented LOSS pseudo-absolute magnitude distribution (see i|3,2H is assumed (Li ct al. 201 lR). 



null hypothesis that the data were drawn from this model 
distribution. In this case, the probabilities for the other 
three distributions are considerably higher, so the model 
does not appear to be inconsistent with the data. 

We also tested the effects of systematically shifting the 
input pAMDs in brightness. If we shift the augmented 
LOSS pAMD 0.3 mag fainter, both the predicted dis- 
tance and absolute magnitude distributions become in- 
compatible with the observed distributions (as above we 
require pipeline detections on 5 or more nights). Shift- 
ing the LOSS pAMD 0.2 mag brighter, the predicted and 
observed distance distributions fall into agreement, but 
the absolute magnitude distributions are not compatible. 
With the SDSS-II pAMD model, a shift of either 0.2 mag 
fainter or brighter can be ruled out. This demonstrates 
that we would be able to detect if the magnitude system 
of the pAMD employed was systematically biased with 
respect to our unfiltered magnitudes by more that ex- 
pected from our checks in £J5J We will use the allowed 
magnitude offsets to estimate our systematic error in 

6. HOST GALAXIES 

iLi et al.l (|2011a[) found that lower luminosity galaxies 



produced more SNe la per unit luminosity than high lu- 
minosity hosts. In order to compare the host luminosity 
distribution of our sample to the SDSS-II and the ex- 
pectations of the LOSS luminosity specific rates coupled 
with a galaxy luminosity function, we have measured 
ugriz photometry for each supern ova host. We make deep 
stacks from the SDSS DR7 data (|Abazaiian et al.ll2009h 
in each of the 5 bands, excluding data obtained while the 
supernovae were active (assumed to be — 20 d< t max < 
100 d in the rest frame), and retain the median using a 
slightly modified version of the Montage packag43- Cuts 
on the sky background and PSF size are used to reject the 
worst images (about a third of th e available data). W e 
identify the host galaxies following lSullivan et al. (2006). 
Briefly, we use the SExtractor shape parameters (C xx , 
C xy , and C yy ) to define an ellipse for each object in the 
field, and we identify the host as the object that includes 
the supernova position inside its ellipse with the small- 
est possible scaling unless the scaling factor, R, is greater 
than 5. If this is the case, we perform forced photom- 

12 http : //montage . ipac . caltech.edu/ 
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Fig. 7. — Distributions predicted from the Monte Carlo simulations (grey) comp ared to the obser ved samp le (green) after correction for 
selection bias. The augmented LOSS pseudo-absolute magnitude distribution (see £|3.2I1 is assumed ULi et al.H2011bl) . 



etry at the location of the supernova in a 3" diameter 
aperture. Otherwise we run SExtractor in dual image 
mode to measure the host galaxy photometry in consis- 
tent aperture sizes on the best matching host. We choose 
to perform the object detection and aperture definition 
on the r-band images. 

A few host galaxies are worthy of special mention: 1) 
The host of SN 2007if is not detected in the co-added 
SDSS data. We inst ead use the g-b a nd m agnitude of 
the host reported by iChildress et al.1 (|2011l ) to scale a 
spectrum of the host galaxy and then measure synthetic 
photometry in the SDSS bands from this. 2) The LOSS 
team assigned the giant galaxy KUG 1259+286 as the 
host of SN 2005ck, but the SDSS reds hift for this galaxy 
is incompatib l e wit h the supernova (|Pugh et al.l 120051 
iLeaman et al.l 1 2 1 lh . There is a small galaxy nearby, 
SDSS J130218.65+282046.8, for which we find a spectro- 
scopic redshift consistent with SN 2005ck; however, this 
galaxy is slightly too far away to meet the host selection 
criteria described above (as would be KUG 1259+286), 
and it is therefore considered hostless. Identified or not, 
the host must be a dwarf, so this should not affect our 



conclusions. 3) SN 2007sp is located a considerable dis- 
tance (R > 10) from a large galaxy at a similar SDSS 
derived redshift. Although t his may yet be the ho st, for 
consistent application of the lSullivan et~a l. (2006) selec- 
tion method, we must consider this target to be host- 
less. 4) There is a galaxy coincident with the location of 
SN 2007kh, but spectroscopic follow-up shows this to be a 
z ~ 0.5 galaxy, which is well in the background. The true 
host of SN 2007kh is therefore a mystery, so we simply 
adopt the photometry of the background source, assume 
the redshift of SN 2007kh, and take this as an upper limit 
on the host luminosity. 5) The location of SN 2007op is 
outside of the SDSS DR7 footprint; however, the host is 
included in D R8, so we adopt th e photometry from the 
DR8 pipeline (jAihara et al.ll201lD . 

In addition, we perform the same analysis on the full 
set of 516 spectroscopically confirmed, probable, or pho- 
to metrically probable z < 0.3 SNe la from the SDSS- 
II iDildav et al.l 120101 The larger number of visits to 
Stripe82 over the course of the SDSS-II results in co- 
added images (again excluding those contaminated by 
supernova light) with limiting magnitudes considerably 
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Fig. 9. — Search efficiency assuming the modified LOSS pseudo- 
absolute magnitude distribution (blue) or the SDSS-II model (red) 
and averaged over all fields and the full survey time-span (including 
un-searched periods). The sudden jumps in efficiency mark the 
minimum allowable distances for various fields (e.g. targets in the 
Virgo fields must be at least 72 Mpc away to meet our selection 
cuts) . 

deeper than the search itself. As a result, even M ~ — 17 
dwarf galaxies can be detected out to the z ~ 0.3 limit of 
the sample. As with the ROTSE-IIIb sample, we manu- 
ally reviewed the selection of each host galaxy and noted 
any special circumstances that might produce mislead- 
ing results (the potential error rate appears similar for 
the ROTSE-IIIb and SDSS-II samples). As needed, we 
adjusted the SExtractor parameter DEBLEND_MINC0NT to 
either ensure large galaxies were not inadvertently split 
into smaller pieces or to attempt to split a likely galaxy 
blend. The latter issue was not always resolved, both for 
the ROTSE-IIIb and SDSS-II hosts.' Of particular note, 
a limited number of SDSS-II hosts appear to be dwarf 
galaxies heavily blended with (if not merging into) larger 
galaxies, such as SDSS-II SN03592. When such blends 
are not separated by SExtractor, we perform forced aper- 



ture photometry at the location of these supernovae to 
better reflect the luminosity of the dwarf hosts. 

Since the host galaxies are distributed over a range of 
rcdshifts, we must perform K-corrections to compare the 
galaxy luminosities in the same rest-frame band passes. 
This is acco mplished using the kcor rect .pro package, 
version 4.2 (jBlanton fc Roweisl [2007h . For comparison 
against the literature, we calculate absolute B- and K- 
band magnitudes for our hosts in the Vega system. These 
are recovered from the best fit galaxy template through 
kcorrect .pro. We compared our derived K-band mag- 
nitudes (based only on the optical SDSS photometry) for 
a fe w low-redshift hosts to detections in the 2MASS sur- 
vey (fSkrutski e et al.ll2006[ ). and found good agreement. 

The four lowest luminosity hosts from the ROTSE-IIIb 
and nearby SDSS-II samples are shown in Figure [TU] and 
the four highest luminosity hosts from each sample are 
shown in Figure 111! In Figure 1121 we compare the host 
galaxies of the ROTSE-IIIb sample to the 20 hosts from 
the S DSS-II sample with redshifts z < 0.09 (pildav et all 
2010). These two samples cover a similar redshift span 
and can thus be compared directly. As can be seen, 
our sample includes relatively few high-luminosity host 
galaxies and, of particular interest, a number of very low 
luminosity hosts. A KS test gives only a 4% chance that 
these ROTSE-IIIb and low redshift SDSS-II samples are 
drawn from the same population. If we assign SDSS-II 
SN03592 to the larger host, this probability drops to just 
1%. The agreement likelihood with the full population 
of galaxies hosting the SNe la from the SDSS-II rate cal- 
culation is about 0.2%. For this larger sample, a small 
fraction of cases that have undetected hosts also have 
limiting magnitudes brighter than detected ROTSE-IIIb 
hosts, so the exact distribution below Mb ~ —17 is un- 
certain. However, this does not affect our key result: 
43% of the ROTSE-IIIb hosts are M B > -18 mag dwarf 
galaxies while only 17% of the full SDSS-II and 2 out of 
20 from the z < 0.09 sample are hosted by such dwarfs 
(including SDSS-II SN03592). 

Figure[T2lals o shows the ho st galaxy distributions from 
the SNLS lAstier et al.l[2006h . As with the SDSS data, 
we measure the host galaxy photometry using SExtractor 
in dual image mode, this time using the i-band as refer- 
ence, and we determine the rest frame absolute magni- 
tudes with kcorrect .pro. We use the MegaPipe me dian 
image stacks provided through CADC (|Gwvni r2008) . To 
check if these measurements are biased by the inclusion 
of images with the supernovae detected in the determi- 
nation of the medians, we compare ou r i-ban d observer 
frame measurements to iSullivan et al.1 ((2010). Surpris- 
ingly, we find that our m easurements are actua lly system- 
atically fainter than the ISullivan et al.1 (|2010| ) values by 
~ 0.3 mag. For example, for the host of SNLS 03Dlar, we 
find vtii = 19. 84 without correcting for 0.05 mag of Galac- 
tic reddening (|Schlegel et al.lH998l ). This agrees with the 
value in the MegaPipe catalog, and it is consistent with 
the SDSS DR8 (19.88 ± .04) but it is fainter than re- 
ported in Sull ivan et al.l (|2010l ) who find nrii — 19.57. 
Aside from this systematic offset, we do not find any 
evidence that our photometry is biased by the super- 
nova light even in the faintest hosts. It is interesting 
that the B-band distributions of the SDSS-II and SNLS 
host galaxies appear to diverge toward the bright end, 
but there is a 12% chance that a random draw from the 
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Fig. 10. — The four lowest luminosi ty hosts from the R OTSE-IIIb SNe la sample (top row) compared to the four faintest hosts from the 
SDSS-II z < 0.15 sample (bottom row; Dilclay ct al. 2010). Each panel shows a 30 X 30kpc close-up centered on the SN position made from 
SPSS images. SN 2007sp is offset from a luminous host at the same redshift, but this separation is more than double that allowed by the 
ISullivan et a l. ( 2010) host selection criteria, so we must consider the host undetected. SDSS-II SN03592 may be hosted by a dwarf galaxy 
merging with a larger host. We were unable to separate these galaxies with SExtractor, so we instead perform forced aperture photometry, 
which may underestimate the brightness of the dwarf. A faint galaxy is seen at the position of SN 2007kh, but spectroscopy indicates this 
is a background source, so we take its measured brightness as an upper limit on the actual host. 



|SN 2008b J 


\ |SN 20Q7sw| 


[s\ 200Sir| 1 


|sN 2006cj| 




1 


# ' 


i 




r|sDSS-II SN17366P 




|~BSPSS-II SN12930|~ 




| 




m 



Fig. 11. — Similar to Figure HOl but for the four highest luminosity hosts in the ROTSE-IIIb sample (top row) and the SDSS-II (bottom 
row). 



14 



Quimby et al. 



1.0 

I 0.8 
a 0.6 

o 

S 0.4 
£ 

0.2 

0.0 
1.0 

a o.s 

q 0.6 

o 

S 0.4 
£ 

0.2 
0.0 



Host Distributions 




ROTSE-IIIb 
SDSS-II z<0.09 
SDSS-I1 (full) 
SNLS 

LOSS (actual) 
- - LOSS (predicted) 
I I I I 



Host Distributions (K-band) 




In this section we calculate the volumetric SNe la rate, 
7£, using: 
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Fig. 12. — Absolute magnitude distributions for the host galax- 
ies of SNe la. A few of the hosts fainter than Mb > — 17 were not 
detected. We employ upper limits for these events, but this does 
not affect our results. The host galaxies of the ROTSE-IIIb super- 
novae (green line) are typically fainter that hosts from the SDSS-II 
search in a similar redshift range (red line), and the ROTSE-IIIb 
sample includes a larger fraction of dwarf hosts than seen or pre- 
dicted through the other surveys shown. 



SDSS-II population would lead to an equal or larger dis- 
placement. As with the SDSS-II, the deep coadds are 
sufficiently deep to detect all SNLS hosts brighter than 

about M B 17, and again the ROTSE-IIIb SNe la 

prefer fainter hosts. 

Although LOSS targets specific galaxies, which biases 
the sample against low luminosity hosts, they adopt the 
IKochanek et al.l ([20011 ) K-band galaxy luminosity func- 
tion and factor in their r ate-size relation to determine the 
volumetric SNe la rate (|Li et al.ll2011af) . We include the 
K-band host galaxy distribution expected by the LOSS 
rate calculation in the lower panel of Figure [12] Clearly, 
this expected distribution is different than the observed 
SDSS-II and SNLS distributions. Specifically, LOSS as- 
sumes that 80% of SNe la in a given volume should come 
from hosts fainter than about K > —23.4 mag, but the 
SDSS-II finds only about 40% of its z < 0.09 hosts are 
fainter than this luminosity. This result holds qualita- 
tively even if the correction for the rate-size relation is 
removed; among the giants, the SDSS-II sample shows 
a preference for the higher luminosity hosts. Moving to 
lower luminosities, the ROTSE-IIIb sample shows a con- 
siderably higher fraction in dwarf galaxies. 

As explained in fj5j our selection process for ROTSE- 
IIIb may have a slight bias against high-luminosity hosts. 
In the next section we calculate the volumetric rates of 
SNe la in dwarfs and giants for comparison against other 
samples. For convenience, we will label galaxies fainter 
than Ms > —18 as "dwarfs," and more luminous galaxies 
as "giants." When discussing the K-band magnitudes of 
the hosts, we will consider Mk > —20 galaxies to be 
dwarfs, since the B — K colors for many of our Mb ~ 
— 18 hosts are around 2 mag. If our sample's apparent 
preference of low luminosity hosts is a result of a selection 
effect, then our SNe la rate in giants should be lower 
than the actual rate while the rate in dwarfs should be 
accurate. We compare our values to published rates as a 
check. 

7. RATES 



n 



eVt 



(2) 



where iV bs is the observed number of events, t is the 
proper time of the survey, V is the co-moving volume 
surveyed, and e is the efficiency factor that gives the 
estimated fraction of events actually discovered to the 
true number of events in the time and volume covered 
(i.e. e < 1). 

We have 18 SNe la appropriate to the augmented LOSS 
pAMD and 16 appropriate for the SDSS-II model that 
were detected on 5 or more nights (see JUand^]). We cal- 
culated the denominator of equation[7]in fJ5]as a sum over 
a series of logarithmically spaced distance bins to account 
for the decline in survey efficiency with distance (sec fig- 
ure [9]) . The proper time for the survey is calculated by 
correcting the observer frame span for time dilation in 
each luminosity distance bin. The survey period was set 
from November 1, 2004 through January 31, 2009, with 
the former date reflecting the month the survey began 
and the latter set to the end of the standard 2008B as- 
tronomical semester. The final term for the survey Vol- 
ume is calculated for each distance bin by integrating the 
co-moving volume element of the Friedmann- Robertson- 
Walker metric with our chosen cosmological parameters 
(a flat, Hq — 71kms~ 1 Mpc -1 , fl m = 0.27 universe) and 
factoring in the fraction of sky covered by our survey. 

Special care must be taken to determine the actual 
survey area since the galaxy clusters of the survey are 
covered by an irregular grid of overlapping fields (see 
Figs. [Hand El). The typical pointing error of ROTSE- 
IIIb is about ±0.05 degrees, so we assume the reference 
templates, which represent the intersection of numerous 
individual images, cover 1.75 x 1.75 degrees each. The 
actual pitch for the main search fields is 1.65 degrees, so 
this is the effective size for fields bordered on all sides. 
To determine the total survey area, we construct an all 
sky register similar to the images in Figures Q] and [U and 
then total all the sky area visited at least once during 
the survey period. Our total sky coverage is 499.0 square 
degrees. 

Adopting the LOSS SNe la pseudo-absolute mag- 
nitude distribution with a supplemental population 
of high luminosity events (assumed to be 1% of 
the total population; see HS.'2\i we find a rate of 
(6.9±i'e) x 10~ 5 SNe IaMpc -3 yr" 1 (statistical error 
only). With the SDSS-II model pAMD, the volumetric 
rate drops to the lower (but statistically consistent) value 
of (4.9t\i) x 10~ 5 SNe IaMpc^yr" 1 ^. These total 
rates and the rates specific to dwarf and giant hosts are 
shown in Figure [T3] and listed in table [31 

For comparison, the SDSS-II found a SNe la rate 
of (2.81±g;||) x 10~ 5 SNe IaMpc~ 3 yr" 1 using their 
z < 0.15 SNe la and assuming the rate to be constant 
out to this limit ([Dildav et alj |2010). Since the SNe la 
rate increases with redshift and even the lowest SDSS- 
II redshift bin extends out to larger distances than the 
ROTSE-IIIb sample, the SDSS-II rate may be slightly in- 
flated with respect to our more nearby rate measurement. 
In any case, the SDSS-II rate is below our 3er lower limit 



ROTSE-IIIb SNe la Census 15 



TABLE 3 
Type Ia Supernova Rates 



Survey Ara.dwarfs Dwarfs Ara, g iants Giants Aia.total Total 

LOSS 

SDSS-II < 0.15) 5 
ROTSE-IIIb (w/ LOSS LF) 6 
ROTSE-IIIb (w/ SDSS-II LF) 4 



0.31 
0.18 
2.31 
1.23 



+0.25 
0.10 
+0.12 
0.08 
1.41 
0.1)2 
1.00 
0.59 



74 
12 
12 



n og+0.43 

2 63+ 34 
4.63+ 178 



3.69 



-1.32 
-1.42 
-1.05 



79 
18 
16 



2.70 
2.81 
6.94 
4.92 



+0.65 
0.41 
+0.35 
0.32 
+2.07 
1.62 
+ 1.58 
1.22 



Note. — Rates in units of 10 SNc Ia Mpc~ 3 yr 1 hi 
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Fig. 13. — Volumetric SNe Ia rates. The right panel shows the 
total volumetric rate fr om ROTSE-IIIb (g reen circles; this work), 
SDSS-II (red diamond; IDildav et all 120101, and the LOSS SNuK 
rate co mbined with a K-band galaxy lumin osity function (blue 
square; Li ct al. 2011a; Kochanck ct al. 2001). The vertical bars 
show the statistical error only. The panels to the left show the con- 
tributions to the total volumetric rate from dwarf and giant hosts. 
The dotted horizontal lines mark the 3<r lower limits for the com- 
plete ROTSE-IIIb samples in each group. The filled green circle 
shows the ROTSE-IIIb rates computed a ssuming the LO SS SNe Ia 
pseudo-absolute magnitude distribution ( Li ct al. 2011b), and the 
open circles show the values derived using th e same SNe Ia dis- 
tribution assumed for the SDSS-II calculation ( Dilday ct ari l2008l , 
20109 an d removing spcctroscopically confirmed SNe Ia in a similar 
manner to the SDSS-II vetting process. 

derived with the augmented LOSS pAMD and may be 
incompatible with the SDSS-II model pAMD as well at 
about the 2a level. We split the samples into dwarf and 
giant hosts and find that while the ROTSE-IIIb rates in 
giants are marginally consistent with the SDSS-II results, 
our rate in dwarfs is significantly higher than the SDSS- 
II with either the augmented LOSS or SDSS-II model 
pAMD. 

To convert the LOSS host luminosity specific SNe Ia 
rate into a volumet r ic rate , we follow similar steps as 
reported in iLi et aLl ()2011al ). but with an important dif- 
ference. We start with the rate-size relation and SNuK 
measurement from LOSS and combine this with the lo- 
cal ga laxy K-band luminosity function Kochanc k et al.l 
(2001). Like the LOSS estimate, the calculation is 
done separately for early- and late-type galaxies, and 
then combined to give the final result. The pub- 
lished LOSS rate is determined by multiplying the lo- 
cal galaxy K-band luminosity densities (jeariy and j\ a t e ) 
by average values of their SNuK, which are determined 
from a weighted average involving the number densi- 
ties from the galaxy luminosity functions. We choose 
to recalculate the LOSS rate through what we believe 



to be a more direct approach. We simply multiply 
the galaxy luminosity function by the rate-size rela- 
tion function and integrate. With our approach, we 
find the volumetric rate implied by the LOSS sample is 
(2.70±g;|f) x 10~ 5 SNe IaMpc" 3 yr" 1 h^. This is con- 
sistent with , if not slightly lower than, the published 
LOSS rate (|Li et al.l 1201 laf ). We compute the statisti- 
cal error through a Monte Carlo calculation including 
the uncertainty in the galaxy luminosity function, which 
gives a somewhat larger confidence interval than reported 
by LOSS. This rate is consistent with the SDSS-II rate, 
which is an odd result given the apparently different 
SNe Ia populations studied (see §3.4[) . Thus, as was 
the case with the SDSS-II rate, the total rates we de- 
rive from our ROTSE-IIIb sample are higher at the 3<r 
level. We also determine the LOSS rate in dwarfs and gi- 
ants by simply changing the integration interval. Again, 
LOSS and SDSS-II derive rates that are consistent with 
each other for dwarfs, giants, and the total host popula- 
tion even though their SNe Ia populations are apparently 
distinct (see further discussion in jjSJ). 

Using cither SNc Ia pscudo- absolute magnitude distri- 
bution, our rate in giants is comparable to SDSS-II and 
LOSS, so there is no evidence that our selection process 
is biased against high luminosity hosts. However, with 
the LOSS distribution our rate in dwarf hosts is signifi- 
cantly higher than either SDSS-II or LOSS, and with the 
SDSS-II model pAMD, our rate is higher at about the 
2(7 level. 

8. CONCLUSIONS 

We have used the background population of super- 
novae discovered by the ROTSE-IIIb telescope over a 
four year period to derive the volumetric rate of SNe Ia. 
Our rate is somewhat sensitive to the SNe Ia pseudo- 
absolute magnitude distribution (pAMD) assumed in de- 
riving our overall search efficiency, and we find that 
recent rates studies have used different, incompatible 
pAMDs. We have calculated our rates using both the 
LOSS pAMD augmented with a small population of high 
luminosity events (1% of the total) and the same pAMD 
model (and similar sample cuts) employed in the SDSS- 
II rate measurement to quantify our distance dependent 
search efficiency. With the augmented LOSS pAMD, we 
find a rate that is more than double published values in a 
similar rcdshift range. This rate is higher than the SDSS- 
II or LOSS rates at 99.9% confidence (ignoring system- 
atics for the moment). Similarly, performing the calcula- 
tion on our "normal-bright" SNe Ia in a manner directly 
comparable to the SDSS-II measurement (including the 
same pAMD model), our rate is still nearly double that 
found by the SDSS-II out to a slightly higher redshift. In 
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this case, our rate is only higher than the SDSS-II value 
at the 2a confidence level, however. 

The LOSS pAMD is utterly incompatible with the 
SDSS-II model, which, in turn, is incompatible with the 
actual volume limited sample from the first year SDSS-II. 
Neither the SDSS-II model nor the LOSS pAMD allow 
for events as bright as SN2007if, so it would seem that 
none of the available distributions account for the full 
population of SNe la. Additionally, the large extinction 
value reported for the photometrically selected SDSS-II 
SN09266 suggests that these distributions may also fail 
to account for the true distribution of host absorptions. 
Down to about -18.8 mag, the LOSS and SDSS-II model 
pAMDs agree, but the LOSS sample shows almost dou- 
ble the relative contribution from fainter events. The 
solid agreement between the LOSS and SDSS-II rates is 
perplexing given that they are apparently sampling com- 
pletely different SNe la populations. 

The systematic errors in our rate measurement are 
dominated by the uncertainty in the SNe la pAMD, in- 
cluding how such distributions derived from (mostly) fil- 
tered photometry apply to our unfiltcrcd data set. Our 
total rate estimate with the SDSS-II model pAMD is 
30% lower than our rate with the augmented LOSS 
pAMD. Part of this is due to the exclusion of SNe 2007if 
and 2008ab from the SDSS-II model rate due to their 
high luminosities (both of these events also have dwarf 
hosts) . The larger fraction of low- luminosity events with 
the LOSS pAMD makes up the additional difference. 
Our search is mainly sensitive to background supernovae 
with pseudo-absolute magnitudes brighter than about 
M R = -18.5 mag, so at least 42% of SNe la are defi- 
nitely missed by our search if the LOSS distribution is 
assumed, but this fraction falls to just 8% if we use the 
SDSS-II model. We use our rate based on the SDSS-II 
model pAMD to constrain the lower bound on the SNe la 
rate and our augmented LOSS rate for the upper bound. 
Adding in a ±0.2 mag systematic offset to account for 
the allowed uncertainty in our unfiltcrcd magnitude sys- 
tem with respected to the (filtered) pAMDs considered, 
our final SNe la rate (including ler statistical errors) at 
a mean redshift of ~z = 0.05 is then between 3.7 and 13.3 
in units of 10~ 5 SNe IaMpc" 3 yr" 1 h% v 

As our sample is spectroscopically complete, we do not 
include uncertainty for unclassified candidates. This can 
be an important concern for other, non-spectroscopically 
complete studies. For example, the SDSS has reportecF^I 
561 spectroscopically confirmed SNe la, but there are at 
least 1070 photometrically probable SNe la uncovered by 
the same survey, and presumably many more transient 
candidates in all including peculiar SNe la, which may 
not pass the ph otometric screening (Dild av et all [2010: 
ISako et al.ll20ll . 

There is also the question of whether all objects tech- 
nically classified as SNe la from their spectra should re- 
ally be grouped in the same physical category. For ex- 
ample, some have argued that some peculiar SNe la, 
such as S N 2005hk, may no t be thermonuclear super- 
novae (cf. iValenti et al.ll2009T >. One event in our sample, 
SN2006ct, bears some resemblance to this subclass, al- 
though we note some differences as well in H2.ll Remov- 
ing this one event would lower our total rates by about 

13 http: //sdssdp62 . fnal.gov/sdsssn/snlist_confirmed_updated 



6%. On the other hand, some events grouped with CC- 
SNe may actually be SNe la in disguise. In particular, 
some Type Iln may be better suited in the SNe la camp 
(|Hamuv et al.l[2^03tlDildav et al.l[20Tl . Over our survey 
period, we discovered three background Type Iln super- 
novae (SNe2006db, 2006tf, and 2008am), but none of 
these have been linked to Type la explosions. 

In using the augmented LOSS pAMD, we have at- 
tempted to measure the total SNe la rate including Hub- 
ble diagram outliers and other peculiar events. This 
should be directly comparable to the rate measured by 
LOSS who include SN 1991T-like, SN 1991bg-like, and 
SN 2002cx-like events in their sample, but this has not 
always been common practice in the literature. Many 
of the previous SNe la rate measurements have explic- 
itly excluded events that deviate from the sub-sample 
of SNe la that are most useful as cosmological probes. 
For ex ample, the recent SN e la rates measured by the 
SNLS (jPerrett et al.l 12012ft are strictly valid only for 
SNe la with light curve widths close to the nominal 
value (0.8 < s < 1.3, but see also lGonzalez-Gaitan et all 
l20ll . SNe la like SN2003fg (a.k.a. SNLS-03D3bb) are 
not discussed. Similarly, the SDSS-II explicitly removed 
their peculiar SNe la, such as the interacting SN 2005gj, 
and the selection cut on the MLCS2k2 light-curve fits 
biases their sample against high luminosity, SN 19 99aa- 
like events (jDildav et aLll20lalKessler et al.l [20091 ) . Be- 
cause of this, the SNLS and SDSS-II rates only reflect 
the frequency for a fraction of the larger SNe la popu- 
lation. As to the high luminosity SNe la excluded by 
SNLS and SDSS-II, however, we note that if these con- 
stituted more than a few percent of the total population 
then they would completely dominate our flux limited 
sample, which is not the case. 

Although the significance of our rates being higher 
than the canonical value is low, it is worth discussing 
how this may have come about. One point of intrigue is 
our excess of SNe la in dwarf galaxies over what would 
be expected from the LOSS or SDSS-II studies. If our 
rate for these hosts is correct, then either the SNe la rate 
size relation increases for the lowest luminosity hosts as 
compared to the range studied by LOSS, or the faint end 
of the galaxy luminosity function used to convert the 
LOSS luminosity specific rate to a volumetric rate may 
be under-counting the actual supply of dwarf galaxies. 

Such correction factors are not required for the SDSS- 
II rate, so the higher frequency of SNe la in low lumi- 
nosity hosts is puzzling. In the full background sam- 
ple, 7 out of 23 of the ROTSE-IIIb hosts are fainter 
than Mb > —16 mag, as are 5 out of the 18 used for 
the rate calculation with the augmented LOSS pAMD. 
These fractions (~ 30%) are significantly larger than the 
1 out of 79 SNe la hosts in the SDSS-II z < 0.15 sam- 
ple that have such low luminosity hosts (hosts brighter 
than M r ~ —15 mag can be detected through our analy- 
sis out to this redshift limit, but only SDSS-II SN13038 is 
fainter) . If this difference is attributed to a selection bias 
against giant hosts in the ROTSE-IIIb sample, then we 
must increase our rate to account for the missing SNe la 
in giants. But to achieve the same low to high-luminosity 
host ratio as the SDSS-II, our already high rates must be 
increased by factors of 4 to 5! 

Another possibility to consider is that the SDSS-II may 
•P'Scmehow be biased against low luminosity hosts, but this 
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is difficult to imagine. A remote possibility is that the 
use of host photo-z's in selecting targets for spectroscopic 
follow-up may boost the likelihood of classifying SNe la 
in giants rather than in dwarfs, which were unlikely to 
have been detected in the reference data available during 
the search. In that case, there would still need to be a 
second bias to reject the SNe la in dwarfs since photo- 
metrically probable yet non-spectroscopica lly confirmed 
SNe la were included in the SDSS-II rates (|Dildav et al.l 

HSU). 

A possible contributor to the relative infrequency of 
dwarf hosts in the SDSS-II sample may be the cut on 
the MLCS2k2 light-curve fits used in selecting photomet- 
rically probably SNe la (jDildav et al.l I2010T ). This cut 
can remove high luminosity SNe la, such as SN 1999aa, 
which tend to be associated with lower luminosity hosts 
(jKessler et al.ll2009l : [SuTlivan et al.ll2006| ). However, only 
one of the five SNe la from our rate sample that fall in 
Mb > —16 mag hosts can be considered highly luminous: 
SN 2007if. Since the SDSS-II does contain a number of 
higher redshift SNe la whose hosts are not detected even 
in our deep coadds, it does not appear that the SDSS- 
II had a selection bias against either "hostless" SNe la 
or SNe la in dwarf hosts, and it is not obvious why such 
hosts are so much more frequent in the ROTSE-IIIb sam- 
ple. 

Another questions is why is there such a larger fraction 
of faint (pseudo absolute magnitude Mr > —18.8) SNe la 
in the LOSS sample as compared to the SDSS-II? Despite 
being targeted, the LOSS galaxy sample spans a range 
of masses and colors, and the full LOSS galaxy sam- 
ple has absolute B and K-band magnitude distributions 
that are consistent with the SDSS-II SNe la hosts in the 
z < 0.09 range (see Fig. [12]). Applying the LOSS rate- 
size rel ation only to ga l axies brighter than Mk < — 20 
in the iKochanek et al.l (|2001l ) galaxy luminosity func- 
tion and selecting only the 48% of SNe la fainter than 
Mr > —18.8 in the LOSS pseudo-absolute magnitude 
distribution, the minimum rate of low-luminosity SNe la 
implied is 0.4 x 10~ 5 SNe IaMpc~ 3 yr _1 In this 

case, the first year SDSS-II z < 0.12 sample should con- 
tain at least 6 SNe la fainter than —18.8 mag. Even 
with the peculiar SN 2005hk and the photometrically se- 
lected SDSS-II SN09266 (with its heavy host absorption 
of Ay ~ 4 mag), there are still only 3 events in the SDSS- 
II z < 0.12 first year sample at such low luminosities. 

It is not clear how the SDSS-II could have missed the 
faint tail of the SNe la distribution (e.g. SN 1991bg- 
like events around Mb ~ —17.2), yet we know of only 
two such spectroscopically confirmed events discovered 
in the three years o f the SDSS-II (2007jh and 2007mm; 
iMosher et all 12012ft and only one possible SN 2002cx- 
like event (SN 2007ie; lOstman et aLl 1201 ID. Addition- 
ally, the Palom ar Transient Factory ( Rau et all 120091 : 
lLaw e t al. 2009J has recently shown that there are per- 
haps even more varieties of peculiar, low-luminosity ob- 
jects that may be spectroscopically classified as SNe la 
(e.g. lSullivan et alJl201ll : lMaguire et al.ll2011h . but these 
are apparently absent from the SDSS-II as well. The 
apparent lack of low-luminosity SNe la could mean that 



there was an unknown selection bias against such objects 
and that the total rate of SNe la is higher than found by 
the SDSS-II. 

The tension between the LOSS and SDSS-II pAMDs 
poses a pote ntially worrisome sys tematic bias for SNe la 
rate studies. iDildav et al.l |2008) argue that the SDSS- 
II was sensitive to even sub-luminous SNe la below 
z < 0.12, so even if their model distribution were flawed, 
this would not have a significant impact on their rate cal- 
culation (but again, the dearth of low-luminosity SNe la 
as compared to the expectations of LOSS is puzzling at 
best and may suggest the SDSS-II has altogether missed 
some SNe la). However, if there is a low-luminosity pop- 
ulation of SNe la that is not accounted for in the SDSS-II 
model distribution, this will lead to an overestimation of 
th e search effic i ency a t the higher redshifts considered 
in IDildav et al.l (|2010| ) , and thus an underestimation of 
the actual rate. In general, rate studies that push their 
samples out to the highest redshifts possible with their 
surveys have an even stronger Malmquist bias and are 
thus more sensitive to systematic errors in their adopted 
pseudo-absolute magnitude distributions. These distri- 
butions may also vary with redshift. It is therefore im- 
portant to securely determine the true SNe la pseudo- 
absolute magnitude distribution from complete, volume 
limited surveys as could potentially be done by the Palo- 
mar Transient Factory at low redshifts and at moderate 
redshifts with the Hyper-SuprimeCam on Subaru. 

Another issue that may lead to bogus rate evolution 
measurements is the changing definition of the SNe la 
sub population measured by different surveys. Low to 
moderate redshift surveys often use l ight curve fits to 
select SNe la (e.g^lDiidav et al.l l2010t INeill et al.l 120061: 
iPerrett et al.l 120121 1 . but the highest redshift studies of- 
ten have far fewer epochs available due to the increased 
follow-up cost, and must t herefore rely more heavily on 
color selection techniques (jGraur et al.l 120 111) . It is not 
clear that these different cuts capture identical sub pop- 
ulations so that differences in the rates between different 
surveys may be taken directly as evidence for evolution 
in the SNe la rate. This question may be resolved by 
performing a rolling search similar to SDSS-II but sensi- 
tive to very high redshift SNe la, as will be possible with 
the Hyper-SuprimeCam Survey. 
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